Experimental determination of a nonclassical Glauber-Sudarshan P function 
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A quantum state is nonclassical if its Glauber-Sudarshan P function fails to be interpreted as a probability 
density. This quantity is often highly singular, so that its reconstruction is a demanding task. Here we present 
the experimental determination of a well-behaved P function showing negativities for a single-photon-added 
thermal state. This is a direct visualization of the original definition of nonclassicality. The method can be 
useful under conditions for which many other signatures of nonclassicality would not persist. 
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Einstein's hypothetical introduction of hght quanta, the 
photons, was the first step toward the consideration of non- 
classical properties of radiation [1]. But what does nonclassi- 
cality mean in a general sense? A radiation field is called non- 
classical when its properties cannot be understood within the 
framework of the classical stochastic theory of electromag- 
netism. For other systems, nonclassicality can be defined ac- 
cordingly. Here we will focus our attention on harmonic quan- 
tum systems, such as radiation fields or quantum-mechanical 
oscillators, for example trapped atoms. 

In this context the coherent states, first considered by 
Schrodinger in the form of wave packets [2], play an impor- 
tant role. They represent those quantum states that are most 
closely related to the classical behavior of an oscillator or an 
electromagnetic wave. For a single radiation mode, the coher- 
ent states I a) are defined as the right-hand eigenstates of the 
non-Hermitian photon annihilation operator a, a\a) — a\a); 
cf. e.g. [3]. A general mixed quantum state p, 

p = j SaP{a)\a){al (1) 

can be characterized by the Glauber-Sudarshan P function [3, 
3]. In this form the quantum statistical averages of normally 
ordered operator functions can be written as 

(: /(a, a^) ■) ^ J d^aP{a)f{a, a*), (2) 

where the normal ordering prescription : f{d,d^) : means that 
all creation operators are to be ordered to the left of all 
annihilation operators a. 

Formally, the resulting expressions (|2]l for expectation val- 
ues are equivalent to classical statistical mean values. How- 
ever, in general, the P function does not exhibit all the proper- 
ties of a classical probability density. It can become negative 
or even highly singular. Within the chosen representation of 



the theory, the failure of the Glauber-Sudarshan P function to 
show the properties of a probability density is taken as the key 
signature of quantumness [5,6]. 

In this Rapid Communication we demonstrate the experi- 
mental determination of a nonclassical P function. Within 
the experimental precision it clearly attains negative values. 
This is a direct demonstration of nonclassicality: the negativ- 
ity of the P function prevents its interpretation as a classical 
probability density. 

Why is it so difficult to demonstrate the nonclassicality di- 
rectly on the basis of this original definition? Let us go back 
to a single photon as postulated by Einstein. Its P function is 

cf. e.g. |7]. Already in this case we get a highly singular dis- 
tribution in terms of derivatives of the S distribution, which 
cannot be interpreted as a classical probability. Due to these 
properties, it is difficult to experimentally determine nonclas- 
sical P functions in general. 

How can one realize nonclassical states whose properties 
can be demonstrated directly in terms of the original defini- 
tion, that the P function fails to be a probability density? This 
question is not trivial: for instance, losses introduced by im- 
perfect experimental efficiencies lead only to rescaling of the 
quadrature variable; cf., e.g., |8J. The P function obtained by 
perfect detection is related to Pri{a), obtained with the quan- 
tum efficiency rj via 

Pia) ^ vPr^i^a). (4) 

Consequently, singularities in the P function are then pre- 
served. Most of the nonclassical states experimentally gen- 
erated so far have highly singular P functions, whose recon- 
struction is impossible. However, one may start with a ther- 
mal state pth with mean photon number n. By photon cre- 
ation one gets a single-photon-added thermal state (SPATS), 
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p = Ma) ptha, where denotes the normalization. Its P 
function is now well behaved, but violates the properties of a 
classical probability density (|9|], 
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P{a)^ — [{\ + n)\a\'-n] 



(5) 



giving rise to the question whether its experimental determina- 
tion could be possible fio']. In the zero-temperature limit, the 
SPATS includes the special case of the one-photon Fock state 
with the highly singular P function given in Eq. ([3]i. In this 
sense the SPATS represents a single photon whose P function 
is regularized by a controlled thermal background. 

Recently, SPATSs could be realized experimentally and 
some of their nonclassical signatures have been verified fTH]. 
Nevertheless, the reconstruction of a nonclassical P function 
remains a challenging problem which goes beyond the stan- 
dard procedures of quantum state reconstruction, for the lat- 
ter, see, e.g., |7]. A successful determination of the P function 
of a SPATS would visualize the basic definition of nonclassi- 
cality for a quantum state that lies at the heart of Einstein's 
hypothesis: a regularized version of a single photon. 

The core of the experimental apparatus used to produce 
SPATSs is an optical parametric amplifier based on a type-I 
/3-barium borate (BBO) crystal pumped by radiation at 393 
nm (see Fig. [TJ. The pump is obtained by second harmonic 
generation in a lithium triborate (LBO) crystal of a mode- 
locked Ti:sapphire laser emitting 1.5 ps pulses with a rep- 
etition rate of 82 MHz. When the parametric amplifier is 
not injected, spontaneous parametric down-conversion takes 
place, generating pairs of photons at the same wavelength as 
the laser source along two directions commonly called the 
signal and idler channels. We perform a conditional prepa- 
ration of the quantum states by placing an on-off photodetec- 
tor (D) after narrow spectral-spatial filters (F) along the idler 
channel 111 ll 11211 . A click of the idler detector prepares the 




Figure 1: Scheme for the conditional excitation of a thermal light 
state (denoted by pin) by a single photon. A click in the on-off detec- 
tor D prepares the photon-added thermal state pout and triggers its 
balanced homodyne detection (BHD). 

signal state, whose quadratures are measured on a pulse-to- 
pulse basis using an ultrafast balanced homodyne detection 



scheme fis'l. After verifying the phase independence of the 
quadrature distributions, the state is then analyzed by acquir- 
ing quadrature values with random local oscillator phases. 
When no fields are present at the inputs of the parametric 
amplifier, conditioned single-photon Fock states are sponta- 
neously generated in the signal channel 1IT2I [Till . On the 
other hand, we have recently shown that the injection of pure 
or mixed states results in the conditional production of their 
single-photon-added versions, always converting the initial 
states into nonclassical ones 

Here we use a pseudo thermal source, obtained by inserting 
a rotating ground glass disk in a portion of the laser beam, for 
injecting the parametric amplifier and producing SPATSs. The 
scattered light forms a random spatial distribution of speck- 
les whose average size is larger than the core diameter of a 
single-mode fiber used to collect it. When the ground glass 
disk rotates, light exits the fiber in a clean collimated spatial 
mode with random amplitude and phase fluctuations, yield- 
ing the photon distribution typical of a thermal source IitIi . 
The product between the SPATS preparation rate and the co- 
herence time of the injected thermal state (a few microsec- 
onds, and depending on the rotation speed of the disk) is kept 
much smaller than 1. This condition assures that each state 
is prepared by adding a single photon to a coherent state hav- 
ing an amplitude and phase which are completely uncorrected 
with respect to those of the previous one. This experimental 
realization of a thermal state directly recalls its P function def- 
inition, i.e., a statistical mixture of coherent states weighted by 
a Gaussian distribution: Pia) = cxp(— |Q;p/n)/(n7r). 

By performing measurements on single-photon Fock states 
and on unconditioned thermal ones, we have estimated an 
overall experimental efficiency of 0.62±0.04. Both the li- 
mited efficiency in the state preparation (wO.92) and in ho- 
modyne detection (wO.67) degrade the expected final state by 
introducing unwanted losses. This does not contaminate the 
obtained P function; cf. Eq. (01). 

Let us now proceed with the reconstruction of the P func- 
tion. Its characteristic function <I'(/3) is related to that of the 
quadrature x[^p) 



$(/3) = (:i)(/3) 



gj|/3|£(V2-arg(/3))\^g|/3|V2^ (g) 



where P>{fi) is the displacement operator Since the measured 
state is independent of phase, we may neglect the arguments 
of and X. The expectation value on the right-hand side repre- 
sents the characteristic function of the observable quadrature. 
It can be estimated from the sample of N measured quadrature 
values {xj\^^^ via (cf. ifT^ ) 
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Inserting Eq. ^ into (|6]l, we get an estimation of $(/3). 
The variance of this quantity can be estimated as 

1 



N 
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(8) 
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The inverse Fourier transform of <I>(/3) yields the P func- 
tion, which for many nonclassical states does not exist as a 
well-behaved function. However, the sampled characteristic 
function converges stochastically toward the theoretical one. 
In our case its Fourier transform is an analytical function. 
For radial symmetry of the state the two-dimensional Fourier 
transform reduces to the Hankel transform il9ll . 

2 f°° 

P{a) = - / 6Jo(26|a|)$(6)d6. (9) 
Ti" Jo 

In our treatment we set the experimental curve to zero for ar- 



The systematic error 

Ap{a) = - / bJo{2b\a\)^b)db (12) 

is estimated with the help of the fitted theoretical function. 

In Fig. |2] we show experimental curves for characteristic 
functions. Curve (a) is in good agreement with the expected 
characteristic function for a SPATS, 

[l-(l + n)|/3|2]e-«l'5|^ (13) 

for the mean thermal photon number n — 1.11 and the global 
quantum efficiency 77 = 0.60. Curve (b) shows the charac- 
teristic function for a mixture of a SPATS and its thermal 
background with weights of 0.81 and 0.19, respectively, for 
n — 3.71 and rj = 0.62. For sampling these functions, we 
have acquired 10^ and 5 x 10^ data points for the curves (a) 
and (b) respectively. We note that both curves are suited to re- 
construct the corresponding P functions by properly choosing 
cutoff values |/3|c of their arguments. 

The reconstructed P function, shown in Fig. [3] is de- 
rived from the experimental characteristic function given in 
Fig. 12 a). Since the measured states are independent of the 
phase, the reconstructed P-representation is phase indepen- 
dent as well. It is clearly seen that the P function attains neg- 
ative values, so that it fails to have the properties of a classical 
probability density. This is a direct proof of the nonclassical- 
ity of the experimentally realized SPATS, based on the origi- 
nal definition of nonclassicality jslSl- 

For a more careful discussion, we also examine a cross sec- 
tion along a radial line, as shown in Fig. SJa). The experi- 
mentally determined curve is drawn with the solid line. Ob- 
viously, it is in good agreement with the theoretical expec- 
tation (dashed curve). The distance between the minimum 
value and the |a| axis is approximately equal to five standard 



guments greater than a cutoff value \(3\c, where the graph be- 
comes small. This limits the disturbing sampling noise on the 
reconstructed function 



_ 2 /■l'^l'= - 

P{a) = - / bJoi2b\a\)^{b)db (10) 
Jo 



to a reasonable level. The corresponding variance has been 
calculated as 



(11) 
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Figure 2: Experimental characteristic functions (solid lines) and best 
fit to theoretical curves (dashed lines): (a) SPATS , with n = 1.11 
and ri = 0.60, (b) mixture of SPATS with 19% of the thermal back- 
ground, with n — 3.71 and 77 = 0.62. The shaded areas show the 
standard deviations. 

deviations, which is not diminished by the systematic error of 
|Ap(a)| < 0.07|P(0)|, obtained by the cutoff \l3\c = 2.8. 
The statistically significant negativity of the P function pre- 
vents it from being interpreted as a classical probability den- 
sity. This provides a clear evidence of nonclassicality per def- 
inition. 

Special nonclassical signatures of SPATSs, which are con- 
sequences of the negativities of the P function, have been ex- 
perimentally demonstrated recently |11]. It is important to 
note that the reconstruction of the P function is just possi- 
ble for sufficiently large thermal photon number n. On the 
contrary, other criteria for nonclassicality, such as negativities 
of the Wigner function, the Klyshko criterion, and the entan- 
glement potential, start to fail for increasing values of n. To 



bb'jQ{2b\a\)Ja{2b'\a\)^{b - b')e''''' dbdb' 



P{af 



4 




■-2. 



0. 

-1. Re(a) 



Figure 3: Experimentally reconstructed P function of a SPATS, as 
obtained from Fig.|2{a). 
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sicality by violating the condition < 1, cf. 12 ill . If 

the condition is satisfied, $(/3) may be integrable and then 
the P function can be obtained to directly verify nonclassi- 
cality. SPATSs belong to this category: for sufficiently high 
n most criteria for nonclassicality (including the lowest-order 
one based on the characteristic function) fail 01 ill , but it is still 
possible to retrieve a negative P function. 

Let us consider how sensitively the negativities of the P 
function depend on the overall efficiency 77. Balanced homo- 
dyne detection measures the "true" state quadratures when 
the efficiency is unity. For imperfect detection (77 < 1) 
one records a convolution of the quadrature distribution with 
Gaussian noise, whose variance increases with decreasing 
77; cf. II22II . In the Wigner function, this increasing noise 
smoothes out its structures and may destroy their negativi- 
ties. As can be seen from Eq. (|4]i, the shape and the relative 
noise level of the reconstructed P function do not depend on 
the efficiency. Hence the negativities of P{a) are in princi- 
ple preserved even for a small efficiency, whereas for other 
phase-space distributions, such as the Wigner function, they 
are quickly lost. 

In conclusion, we have reconstructed the Glauber- 
Sudarshan P function of an experimentally prepared single- 
photon-added thermal state. We obtain a well-behaved func- 
tion with statistically significant negativities, so that it fails 
to show the properties of a classical probability density. This 
is a direct demonstration of nonclassicality according to its 
original definition. The approach works well, just when many 
other methods of demonstrating nonclassicality fail. 

This work was partially supported by Ente Cassa di 
Risparmio di Firenze and CNR, under the RSTL initiative. 
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Figure 4: The P functions (solid lines) in parts (a) and (b) are ob- 
tained from the experimental characteristic functions in Figs. |2la) 
and|2fb) respectively. They are compared with the corresponding 
theoretical fits (dashed curves). The standard deviations (light shaded 
areas) and the systematic errors (dark shaded areas) are also given. 



show the power of the reconstruction of the P function under 
such conditions, we have demonstrated its use at the limits: 
for a SPATS with n = 3.71, which is additionally contam- 
inated with a 19% admixture of the corresponding thermal 
background. By using a cutoff \(3\c ~ 1.9, we still obtain 
a P function being negative within one standard deviation. 
Other nonclassical effects, as discussed above, do not survive 
for this state. 

Criteria for the characteristic functions are known, which 
are equivalent to the negativity of the P function i20ll . For 
many states the characteristic function displays their nonclas- 
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